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N. M. Reddy» 

Langley Research Center 


SUMMARY 

Convective heat-transfer measurements were made on the conical 
portion of spherically blunted cones (30° and 1*0° half angle) in the 
Langley expansion tube. The test gases used were helium and air; flow 
velocities were about 6.8 km/sec for helium and about 5«1 km/sec for air. 
The measured heating rates were compared with csdculated results using 
a viscous shock -layer computer code. For air, various techniques to 
determine flow velocity yielded nearly identical results, but for helium 
the flow velocity varied by as much as 8 percent depending on which 
technique was used. The measured heating rates were in satisfactory 
agreement with calculation for helium test gas if the higher flow 
velocity was used, but for air and for helium assuming the lower flow 
velocity, the measurements were significantly greater than theory and the 
discrepancy increased with increasing distance along the cone. 


INTRODUCTION 

Determination of windward side heating rates during entry of the 
Space Shuttle Orbiter is very difficult. Analyticsil methods for 
calculating flow fields over complex configurations are as yet in a 
developmental stage, are costly, cumbersome, and limited in scope of 
variables. The extent of experimental data for such configurations is 
also limited, particularly for high velocity flow conditions. One 
approximation which has been made to permit calculation of the windward 
heating distribution on a vehicle at high angle of attack is to assume 
the flow along the windward centerlino is similar to that of a cone 
whose half-angle is equal to the vehicle's local slope relative to the 
velocity vector. Several methods ore available to the engineer to 
calculate the flow and heat -transfer rate for^ the equivalent blunted 
cone at r.ci*o angle of attack. The critical heating rates on windward 
surfaces of the space shuttle are expected to occur at altitudes between 
6l and 76 kilometers (200,000 and 250,000 ft,). In an effort to 
simulate the shuttle flow field, total flow-field calculations for blunt 
spherical nose cones (half-angle 30° and 1*0°) have been presented in 
reference 1. This analysis consisted of defining the complete flow field 
around blxinted cones for assumed equilibrium and reacting gas cases. The 
axial length of the blunt cones was assumed to be 30.5 meters (lOO ft.). 
This analysis (ref. l) was mainly concerned with delineoting in detail 
the real-gas effects in the boundary layer and their effect on the surface 
skin friction and heat transfer. 


*NRC-Senior Research Associate, on leave from Indian Institute of Science, 
Bangalore, India. 



The simulation capabilities of various Langley facilities have been 
also assessed in reference 1. The facilities considered are the Mach 8 
variable density wind tunnel, Mach 20 helium tiuinel, Mach 20 nitrogen 
tunnel, and the Mach 6 Freon (CFi^) tunnel. However, no simulation 
capability of the expansion tube facility is assessed in reference 1. 

The results showed that the best simulation of the Invlscid flow field 
(shock standoff distance, edge Mach number) was provided by the CFi* 
facility. The laminar heating simulation provided by all facilities was 
eqvd. valent. It was concluded in reference 1 that natural boundary -layer 
transition is unlikely except in the variable density wind tunnel, 
especially at higher stagnation pressure level. Based on the results 
presented, the authors of reference 1 have concluded that investigations 
of two types would be required to verify the analyticeJ. results, namely 
(a) basic investigations of boundary -rlayer profiles combined with skin 
friction and heat -transfer measurements on relatively simple cylindrical 
shapes, and (b) aerodynamic heating tests on blunt-cone configurations of 
the type analyzed in reference 1. 

The present investigation measured the heating-rate distributions 
over blunt cones of 30° and 1*0° half-angle in the Langley expansion tube 
facility with helium as well as air as test gases. The purpose of the 
present study is to present hypersonic cone heating data at both perfect* 
gas and real-gas conditions. The measured heating rates have been 
compared with computed values (ref. 2). An attempt has been made o 
racusure heating rates within an accuracy of 1* to 5 percent by iising two 
different and independent methods to deduce heating rates. Based on the 
results of the present investigation, some conclusions are drawn regarding 
the accuracy of heating -rate measurements and their comparison with the 
computed values. 


SYMBOLS 

speed of sound in driven tube 

c coefficient of specific heat at constant pressure for quartz 

^ substrate 

C capacitance in the analog network 

D amplification factor for analog circuit 

E voltage across the thin film and also across the analog circuit 

AE voltage change across the thin film 

initial voltage across the thin film 
Ejj battery voltage 

1 current through the gage 
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initial current through the gage 
current through the analog network 

coefficient of thermal conductivity’ of quartz s\fl>strate 
incident shock wave in the i) termediate section 
incident shock wave in the acceleration section 
free-stream Mach number 

flow Mach number behind the incident shock wave M„ . 
free-stream pressure 
pressure behind shock wave 

initial pressure in the intermediate section or driven tube in 
the calibration shodc tube 

Prandtl number 

initial pressure in the acceleration section 
heat-transfer rate (watts/m ) 
stagnation-point heat-transfer rate 
expansion tube radius 
body nose radixis 

resistance in the analog network (also gage resistance in 

appendix C) 

defined in appendix C 

initial resistance of gage 
free-stream Reynolds nundier /meter 
distance along the body 
time 

free-stream temperature 
temperature of quartz substrate 



Uo 


u. 


u 


s ,10 


u, 


s,l 

27-31 


U 


average 


velocity at the boundary-layer edge at the stagnation point 

flow velocity behind shock wavi? M - 

s ,1 

flow velocity in region (fig. l) 

shock velocity in the acceleration section 

shock velocity in the intermediate section 

shock velocity determined between station numbers 27 and 31 

shock velocity at model location from average curve fit 


V voltage across the first resistance in the analog network 

X distance from tube centerline normal to flow axis 

z distance from tube exit 

e 

a thin-filra temperature coefficient of resistance 

l/P 

0 defined as 3 = (kpC ) ' 

P 

p density of quartz substrate 

e,rt density ratio across the shock M 

10 8,10 

p density at the boundary -layer edge at the stagnation point 

p viscosity at boimdary-layer edge 

G 

0 cone semi vertex angle 

X dummy variable (eq. (l)) 


FACILITY AND TEST CONDITIONS 


Expansion Tube Facility 

Several papers (refs. 3 to 9) have been published concerning the 
description and performance estimation of the expansion tube facility 
at the Langley Research Center. It consists of a driver section 
(2.1*h-m long and 0.l65-m-dia.), an intermediate section (7.L9-m long 
and 0.152-m dia.) and an acceleration section (lU.12-m long and 0.152-m 
dia.). Ihe driver section is pressurized at ambient temperature with 
a gas having a high speed of sound. For the present tests, hydrogen 
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was \i3ed as the driver gas. The Intermediate section Is evacmted 
and filled with the desired test gas with pressure at ambient 
temperature. The acceleration section was also evacuated along with 
the dunp tank to a very lew pressure and filled with the same gas as 
the selected test gas but to a much lower presexire level pj^Q at room 
temperatiu*e. The pressure difference between the intermediate and 
accelotation sections was maintained by a thin secondary diaphragm made 
of polyester film ranging in thickness from 3.175 to 25.^* pm. The main 
diaphragm between the driver and the intermediate sections was made of 
stainless steel. The test models were positioned 0.038 m dewnstream from 
the acceleration tube exit. Flow throu^ the acceleration section 
exhausts into a dunp tank; hence, models are tested in an open Jet. A 
detailed description of the basic components and auxUlary equipment of 
the Langley expansion tube is presented in references 9 and 10. 

The operating sequence for the ejqjansion tube is shown schematically 
in figure 1. The operating sequence starts with riq>ture of the high- 
pressure primary diaphragm. This creates an incident shock wave which 
travels down the intermediate section \mtil it encounters the secondary 
diaphragm. Ihe low^ressure secondary diaphragm will be ruptured due to 
impact of incident shock wave. This creates a secondary incident shock 
wave which propagates into the low-pressure acceleration gas while an 
upstream expansion wave moves into the test gas. In passing through 
this upstream expansion wave, which is being washed downstream since the 
shock-heated gas is siq>ersonic, the test gets tindergoes an isentropic 
unsteady expansion resulting in an increase in the flow wlocity and Mach 
niimber. 


Free -Stream Flow Calibration 

Flow uniformity .- Extensive flow calibration in the expansion tube 
facility for the case of helium test gas has been recently completed and 
reported in reference 9. Pitot pressure profiles at various distances 
downstream of acceleration tube exit in the case of helium test gas 
(Pi ® 3^^7 N/m^; p^^* = l6.0 N/m^, t = 200 ysecs) are shown in figure 2 
(taken from ref. 97. From this study it was found that for the initial 
conditions used in the present tests, the pitot sxurvey shows uniform 
flow region of about 7 to 8 cm in diameter between 1 and 20 cm downstream 
of the acceleration tube exit. Calibrations for air test gas have 
shown a similar test core size. 

Free-stream flow quantities .- Wie free-stream flow quantities were 
computed from the real-gas computer program of reference 11. At least 
three flow quantities (either predicted or measured) have to be known 
to obtain other flow quantities from the computer program of reference 11. 
Calculated flow quantities based on test section flow measurements arc 
considered to be superior to those based on initial state variables 
because of deviations from theoretical prediction. . The usual procedure 
would be to use the measured wall static pressure euid shock velocity near 
the acceleration tube exit in each run and the measxu'ed pitot pressure 
in the test section in a repeated run as the three required flow 
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quantities to compute the other free-stream flow quantities from the 
computer program of reference 11. In the expansion tube facility the 
repeatability of runs as ascertained from the measured shock velocity 
is good in almost all cases with air and helium as test gases. But 
the inherent nonideal effects (viscous, real-gas and shock-reflection 
phenomena at the secondary diaphragm) might affect the flow quantities 
in the test section significantly in each run. In this sense, the 
repeatability as ascertained from the measured shock velocity may not 
be sufficient. Therefore, in the present investigation it was decided 
to measure simultaneously the pitot pressui-e in addition to wall static 
pressiire and shock velocity in each run. In the present investigation 
the pitot pressure was measured by installing a small pressure transducer 
(2.5-mm dia. ) in the blunt nose of each of the heat-transfer cone models. 
Since the transducer size was very small, the effects on flow over the 
model due to curvature change near the nose region were assumed to be 
negligible . 

Free-stream static pressure .- The wall static pressure, which was 
measured at a station 0.026 m upstream of the acceleration tube exit, 
by using a commercial piezoelectric pressure transducer, was assumed 
equal to the free-stream static pressure (ref. 9). 'The transducer was 
connected to a charge amplifier, -he output of which was displayed on an 
oscilloscope. 

Pitot pressure .- The pitot pressure was measured by using a 
piezoresistive pressure transducer with a fully active semiconductor 
Wheatstone bridge which has a ncminal sensitivity of aboxit 0.0725 mv/ 
kN/ra^ (0.5 mv/psi). These gages were calibrated before each run by 
using a step pressure input and the pitot pressure variation during 
each run was recorded and stored at 2.!>-)iS('c intervals in a waveform 
recorder. 

Free-stream velocity .- The free-stream velocity U5 in region 5 
(fig. 1) is ideally equal to the acceleration gas/test-gas interface 
velocity. This velocity is normally not measured but deduced from the 
measured shock velocity ^3,10* shock tube theory shows that 

the separation distance between the shock and interface increases 
linearly with distance fix>m the diaphragm station. However, the shock 
tube wall boundary layer caiases the incident shock to decelerate and the 
Interface to accelerate and thus approach the same value. It is shown 
in reference 9 that for helium test gas, for values of pj^Q less than 
50 N/m^, the interface velocity is essentially equal to the incident 

shock velocity This assumption was made in the present 

investigation. In the case of air test gas also the same assumption 
(U5 - Ug 10^ made since the p^q value used was 6.7 N/ra^. 

The dominant fact here is that the maximum separation distance between 
shook and Interface has been approached for the p^^Q range understudy. 
From the preceding discussion it is obvious that an accurate measurement 
of shock velocity near the acceleration tube exit is essential for an 
accurate prediction of test section flow velocity. 
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The shock velocity Ug i*' "the acceleration section is measured 
routinely in the Langley expansion tube raeility by time-of-arrival data 
gages mounted in the tube wall. More details are foimd in references 
10 and 12. The time-of-arrivtil data were obtained by counter readings 
which were triggered by signals obtained from wall-mounted thin -film heat- 
transfer gages along the length of the acceleration section. Typical 
variations of shock velocity Us,10 helium and air test gases are 
shown in figurts 3 and U, respectively. A straight-line fit to the 
mesisured data is also shown in figures 3 and 1*. The velocity at the tiibe 
exit given by the strai^t-line fit is given in table I for all the rxms. 
Con^jarisons of figures 3 and h show that helitim exhibited more scatter 
in the velocity measurements » because of the lower signal levels. Also, 
some counters in several runs gave completely erroneous readings, 
especially in helium, and these readings had to be omitted in fitting 
a strai^t line to the data. 

Because of the importance of an accxurate velocity determination at 
the test section, an examination was made of the method for determining 
velocity. The shock velocity between the last two stations (station 
numbers 27 and 31 which are 1.53^ m apart) was deduced by using corre- 
sponding counter readings as well as the time interval measured from 
the oscilloscope film of heat-transfer gage outputs. These values are 
also shown in table I. The wall pressure gage film data between these 
two stations were also available and velocity deduced from tliese data is 
also shown in table I. The time sweep of all scopes was calibrated with 
a time mark generator and corrections for sweep rate were applied 
wherever necessary. It is apparent from table I that the shock velocities 
obtained by different methods for any given run do not agree with each 
other. 

In the case of air test gas , the shock velocities deduced from the 
counter reading between stations 27 and 31 are believed to be siispect 
since the output from the heat -transfer gage at station 27 went 
negative immediately after the shock passage before it recovered back 
to its normal behavior and this might have affected the counter reading. 

The discrepancy noticed between velocities obtained from heat transfer 
and wall pressure film data in table I is somewhat surprising since the 
method used is the same. However, it was observed that the heat-transfer 
gage's re.sponse at stations and 31 was much slower than the pros.surc 
gage's response at the same stations, which is perhaps due to thick coating 
on the heat-transfer gages. For this reason, the velocity obtained from 
the pressure film data between stations 27 and 31 nay be of superior 
value. It is interesting to note that in the ease of air, the 
differences between the velocities obtained from the straight-line fit to 
the data over the entire tiibe length and the pressvire film data are very 
small. 

In the ease of heliunw the differences between velocities obtained 
from different methods are much more than those obtained in air. The 
output from the heat-transfer gage at station 27 in the case of helium 
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did not go negative, unlike the case of air, probably d\ie to negligible 
ionization. But the response vas too slow and the signal level was 
very small which might have affected the counters significantly. 

Because of these discrepancies, the theoretical heating predictions 
were con 5 >uted for both the velocity deduced from the curve fit 
and that obtained from the pressure data {U27~31^ helium ^ 

nms. For air runs, the difference was not significant so only 
was used to confute heating rates. Flow^state calcvilatlons were 
based on 

FABRICATION AND CALIBRATION OF THIN- FILM GAGES 

Thin-film heat-transfer gages were fabricated by sputtering palladium 
on quartz substrates. The quartz substrates were fabricated with proper 
radius of curvatxire to fit the cone surface. The palladl\im film was 
sputtered over the entire top surface of the substrate and then the thin 
film in a serpentine pattern (fig. 5) was etched by using a photoetching 
process. The ends of the thin film «,re Joined by a thick film in the 
sides and wires were soldered to the thick film. During development of 
this technique it was found that a thin coat of chromium (sputtering time 
30 secs) on the substrate before palladium application inproved the thin- 
film adherence to svibstrate enormously. The sputtering time for palladium 
film was 5 to 6 minutes. A coat of silicon dioxide (sputtering time 7 
to 8 minutes) was sputtered on top of the palladi\am film to avoid possible 
shorting across the thin film due to ionization in the test gas. About 
lU to l6 gages at nearly 3 mm apart were fabricated on each substrate 
(fig. 5). The quality control was such that all the films were identical 
in their geometrical characteristics; the resistance of each film was 
within a few percent of each other. The resistance value varied between 
75 and 100 0 for different batches. 

Two types of calibrations are required for any thin film to be 
used as a quantitative heat-transfer measuring device* namely, the 
temperature coefficient of resistance a and the substrate property 
S = (kpcp)l/2^ Xhe temperature coefficient of resistance a(ohm/ohm-k) 
was determined for each film at four different tenperatures by heating 
the substrates in a controlled furnace and the a deduced by a 
computer program. Typical data for a substrate given in table IT, 

This method is believed to give o with an accuracy of 2 to 3 percent. 

The standard deviation in a between different gages is within 0.6 
percent as shown in table II. This confirms the high degree of quality 
control maintained during gage fabrication. 

The vtilue of 3 for the quartz substrate cannot be obtained from 
the available bulk physical properties of quartz since uncertainties 
that exist in the physical properties evaluation are unacceptable in 
the present type of investigation. Furthermore, it is well established 
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that the physical behavior of quartz substrate InDSdiate^ behind the 
thin film will be different than the bxilk substrate due to the presence 
of thin film which is actually fired on to the substrate surface. 

Iberefore, several techniques (refs. 13**15) are available to determine 
the value of B after the thin film has been fabricated on to the sub- 
strate surface. Among the various techniques available* the most accurate 
and direct method would be to \ise a small shock t\d>e (ref. 13) and 
measure the stagnation-point heat transfer at low shock Mach nuni>er8 
(1.3 i <. 2). Then using the available theories for stagnation point 
heat transfer (refs. l6 and 17) » the value of 0 can be deduced vezy 
easily. The details of the present calibration method are given in 
appendix A. For quartz substrates* it has been well established that a 
valu<» of 6 = 1512.5 W-sec^/2/5,2„jt used. Prom the present calibration 

data (given in appendix A), the value of 0 wotild be 1532.4 W-sec^/2/m2_jt vhlch 
is about 1.3 percent more than the value assumed by other research workers. 

Since the standard deviation between repeated runs in the present 
calibration (Appendix A) is about 2.5 percent* it was decided to use the 
value of B = 1512.5 W-sec^/2/jn2«j{ jjj present investigation. 


BLUNT NOSE-CONE NOBELS 

Two 4o^ and 30^ semi vertex angle ccxie models with a spherical blunt 
nose (6.35-mm radius) were fabricated out of stainless steel. Ihe 
geometrical details of these models are shown in figure 5* Quartz 
substrates of 6.35 mm x 6.35 mm cross section were fitted in a slot cut 
along a single ray in the conical porticsi of the models as shown in 
figure 5. The top surfaces of the quartz substrates were optically 
ground and polished to match the vaxying curvature along the cone surface. 
About l4 to 16 thin-film gages at intervals of 3.2 mm were fabricated on the 
quartz substrate as shown in figure 5* Photographs of the 46^ cone model 
with quartz substrate and thin films before and after the run are shown 
in figures 6 and 7* respectively. Ihe models were tested in the expansion' 
tube facility at zero and positive angles of attadc with air as well 
as helium test gas. The acceleration gas was the same as the test gas 
in all the runs. Ihe initial conditiois used in the intermediate and 
acceleration sections* as well as the different angles of attack used 
in all the runs* are ^ven in table III. All the thin films were damaged 
by the fragments of secondary diaphragm during each run. Small white 
fragments sticking to the gage substrate are visible in figure. 7* The 
performance of the thin film was not affected during the first 1/2 to 1 
millisec time duration since the secondary diaphragm fragments arrived 
at the model after about 1 millisec ^ich was confirmed by monitoring 
the gage response for a few millisecs dxiring the run. 


HEAT-TRANSFER MEASURING TECHNIQUES 

The heating Hrate distributions over cone models were measured 
simultaneously by using two methods which are coaq;>letely independent and 
different from each other. Some ing>ortant details of these methods are 
given below (more details in appendices B and C). 
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Heating Rates Obtained by Numerical Integration 

The surface-tenqperature time history of the gage s\ibstrate materlol 
(quartz) was measured by using an electric circuit to determine gage 
resistance (details given in appendix C) and the heating rate vas deduced 
by numerical integration of surface temperature history. The relation- 
ship between the heat transfer per unit area per unit time q(t) and an 
arbitrary surface -temperature history AT(t) with AT(0) « 0, is 
(ref. l8): 



1/2 

B s ^He quartz substrate and X is a dummy variable. 

Equetion (l) has been modified for easy integration in reference 19 
and that expression vas used for the numerical integration. The surface- 
temperatiure time histories, measured by using the electric circuit 
(appendix C), were recorded and stored in the waveform recorder and then 
nximerically Integrated in a desktop cedculator. More details are given in 
a later section. 

Heating Rates Obtedned by Electric Analog Network 

In the derivation of the expression given in equation (l), several 
approximations are made (for details, ref. l8 may be referred). The 
numerical integration of equation (l) hao to be performed in each case 
and this requires accurate surface-teinperature data at very close intervals 
of time. In the early sixties, Skinner (ref. 20) a::d Meyer (ref. 21) 
independently realized that electrical analogs of a heat-conducting 
slab could be constructed which, if driven by a voltage proportional 
to surface temperature, would yield a si^al proporzimal to the heat- 
transfer rate. This electric analog not oily avoids the tedious numerical 
integration procedure but also avoids the assumptions made in arriving 
at the expression for heat-transfer rate given in equation (l). 

Purthermore, numerical integration of surface-teo^rature history can 
give erroneous values of q(t) if very accurate values of measured 
AT(t) are not \ised at very close intervals of time. The basic principles 
of electric analog network and the actual circuit diagram used are 
given In appendix B. The relation between heating rate, the analog 
parameters, and substrate properties is given as: 


q(t) 


2V^ 




( 2 ) 
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V^(t) is neastured as an output ffom the analog and all others are known 
cmstants defined in appendix B. Hence* q.(t) is directly proportional 
to V^(;}. In the present investigation 10 channels of electric analog 
networks were designed and bxillti and these were iised to neasxire 
directly the heating rates. The outputs frost the analog channels vetre 
recorded and stored in the waveform recorder and the heating rates 
were directly printed out and plots were made using a desktop calculator. 

Data Recording and Processing 

Heat-transfer measurements over the cone surface were made with cone 
models mounted in the expansion tube test section (model nose located at 
38 mm from the acceleration tube exit) with helium and air as the test 
gases. The initial flow conditions that were \ised are tabulated in 
table III. 

For some of the runs the heat-transfer rates for all the gages on 
the cone surface were deduced from the me8S\ured swface teoiperaturcs by 
numerical interaticn. The surface temperatures were recorded using an 
electric circuit described in appendix C. Since this circuit had no 
preamplifier, the signal was amplified by tising a differential amplifier. 
Proper ground connections to isolate ground loops were essential to avoid 
noise in the surface-temperature histories. The amplified signals were 
recorded on oscilloscopes as veil as in the waveform recorder at 2.5- 
psec intervals in all the runs using 15 out of 17 available channels. 
After each run, the stored data were transferred to a desktop calculator 
and printout of the data as well as plots was obtained. The data were 
also stored in a floppy disk and \ised for numerical integraticsi to 
obtain the heat-transfer rate. Chie of the 15 channels in the waveform 
recorder was used to record the pitot pressure measxired by a pressure 
transducer mounted at the blunt nose, as shown in figure 5. By using 
the proper calibratiMJ factor for the pressure gage, the printout from 
the calculator as well as the plots was obtained in actual pressure 
reading (either in pal or N/m^;. 

For some runs, the heat-transfer rates were directly measured by 
using electric analog networks (details in appendix B). Since the circuit 
used in this case contained an anqplifier in addition to the analog 
network, the signal output from tl'e analog was directly recorded on 
oscilloscopes as well as in the waveform recorder, without any further 
amplification. In this case also, the data were recorded at 2.5-psec 
intervals and processed in the calculator, since the nunber of analog 
network channels was «ily 10; the remaining U to 5 gages on the model 
were connected through the electric circuit (appendix C) and heating 
rates were obtained from the recorded surface tenqperaturea by numerical 
integration. However, for the data obtained from walog networks, 
by using a proper constant (obtained from eq. {2))t the heating rates 
were directly printed out and plotted. Typical plots of surface 
temperatures and heat-transfer rates are shown in figure 6. 





Data Reduction Proced\ire 


Comparison of heating rates measured by analog network and deduced 
from numerical integration .- Tvo types of tests were made to cooipare 
vhe two different methods that have heen tised in the present analysis 
for heat-transfer measurements. Urst* the surface -temperature history 
obtained due to incident shock heating in air in the acceleration section 
was recorded by installing a thin-fllm gage inside the shock -tube wall. 

The surface-temperature history was recorded on a scope and the sane 

signal was simultaneously fed into the analog network and the output from 

the analog was measured in the wave-form recorder. The sm fact— temperatiire 

history was measured from the oscilloscope film using a digitizer and the 

heating rates obtained by this procedure and those directly meastired from t 

the analog are compared in fisairc y. The significant and important • 

observation :.s that the numerical integration technique does not seem to 

simulate the fast response of the heating rate immediately behind the shock 

wave whereas the analog circuit shows a much better response in this region. 

About 35 to 1*0 ysecs after the shorV pnoaage.the two methods agree within 
8 percent. ; 

The heating rates were also deduced \rsing the preceding procedure 
in the case of stapiation-polnt heating-rate measurement to a spherical 
nose of 25<^*inm radius mounted in the expansion tube test section. 5be 
test gas was helium and the other conditions correspond to run number EC 
1500 given in table III. The surface-temperatxore history was recorded 
in the recorder at 2.5-Wsec intervals. The numerical integration was 
performed by using data at two time intervals, namely 2.5 and 5«0 ysecs. 

The heating rates variation with time is shown in figure 10. The 
differences in heating rate by using the two time intervals is negligible. 

However, the computer time required to perform the integration over a 
200-ysec time period was reduced by 50 percent for 5-Wsec time intervals. 

Tlje heating rate obtained from the analog network by feeding the same 
surface temperature was also recorded in the recorder at 2.5 ysec 
intervals ajad is compared in figure 10 with the values obtained by 
numerical integration. In this case also, the numerical integration 
does not seem to simulate the fast response of the heating rate immediately 
behind the shock wave. The two methods agree well (differences within 
3 to U percent, except for a couple of neak values) after an elapsed 
time period of about 90 to 95 ysecs. This fact was taken into account 
In the dats-reduction procedure adopted for the heating rates over 
cone models. 

Heating-rate distributions over cone models .- Tbe measured heating 
rates over the cone models ideally should remain constant with time. 

Tbis ideal situation requires that the frec-stream flow quantities 
are constant with time during the testing period. However, due to several 
nonideal effects that exist in the operation of the expansion tube 
facility, the free-stream flow quantities do exhibl’' some variation 
with time, especially in the case of helium test gas. Therefore, a 
time-averaging procedure was adopted to obtain the measured heating 
rates over the cone models. 
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As discussed earlier, the three required fyee-stream quantities, 
namely, wall-static pressure, pitot pressure, and the flow velocity 
are measured in each run and Ihese are used as inputs to compute the 
theoretical heating-rate distribution. Therefore, it is important 
to obtain a consistent, time-averaged set of measured free-stream flow 
conditions. For every run, the variation of pitot pressure with time 
was plotted from the recorded data. The wsdl static pressure was read 
from the oscilloscope film in a digitizer and plots were made. The 
measured heating rates for two gage locations (genereilly iu the 
beginning and toward the end of the quartz substrate) are also plotted, 
laical plots are shown in figure 11. All three plots are matched with 
time t = 0 and then a common time interval was located in all three 
plots, during which period the flow quantities variation with time is 
minimum, as indicated in figure 11. In selecting this time interval, care 
is taken such that the initial period from t = 0 to the beginning of 
the time interval is at least 80 to 90 ysecs, especially in the case of 
numerical integration, in order to avoid the period of sluggish response 
as already discussed. One exception to this general rule has been made 
in the case of helium, due to nonavailability of better time intervals. 
Tliis specific case has been noted in the discussion of results. After 
selecting the proper time interval, the pitot and static pressures were 
time-averaged. These values, along with the corresponding meaisured flow 
velocity for all the runs, are listed in table IV. The other flow 
quantities obtained from the computer program of reference 11, are listed 
in table V for all the runs. The heating rates were also time-averaged 
during the same time interval and the resulting heating-rate distribu- 
tion along the length of the cone was obtained for sdl the runs. 


THEORETICAL HEATING-RATE DISTRIBUTIONS OVER BLUNTED CONES 

Theoretical heating-rate distributions have been obtained by using 
the method developed by Kumar and Graves (ref. 2) for the viscous flow 
over blunted cones at zero and small angles of attack. This method uses 
time-dependent viscous-shock-layer-type equations to describe the flow 
field between the bodj' and the shock wave. A time-dependent finite- 
difference technique is lised to solve the equations in the planes of 
symmetry of the flow fiell. Complete details of this method are given 
in reference 2. 

Although the above analysis provides solutions at small angles of 
attack also, it is shown by Kmar et al (ref. 22) that a spherically 
blunted cone at small angles of attack can be replaced by eqviivalent 
cones at zero angle of attack for the purpose of calculating the 
shock shapes, surface pressure, skin friction, and heating rates. As 
an example, a 50° sphere cone at 5° angle of attack can be replaced 
by a 55° sphere cone at zero angle of attack on the windweurd side and 
a 1»5° sphere cone at zero angle of attack on the leeward side. This 
approach is used in the present investigation to obtain theoretical resvilts 
for comparison with the measured heating rates at angles of attack. 
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In the present investigation, heating rates over the 30° and 
sphere cones are computed iinder perfect-gas approximation for helium 
test gas and under chemical-eq\ii librium approximation for edr-test 
gas. The viscosity for helium is obtained by the expression 


U = 3.92 X 10"^T°*^ N-sec/m^ 


( 3 ) 


where T is in °K. 

The equilibriiuD meiss fractions of various species of air are computed 
by using the free-energy-minimization technique and their enthedpies are 
obtsiined by using the expressions given in reference 23. The flow 
properties predicted by this method agree within 1 percent from those 
predicted by the computer code of reference 11 at the stagnation point. 

Tlie transport properties for equilibrium air are taken from reference 2k. 

The radiative-heating terms were also included in the governing 
equations. The radiative heat flux model used was that given in reference 
25. In all the coses computed, it was found that the radiative heating 
rate was almost negligible. 

The heating rates for the present experimental conditions were 
kindly computed by Dr. AJay Kumar, Research Assistant Professor, Dept, 
of Mechanical Engineering and Mechanics, Old Dominion University, 

Norfolk, Virginia. 


RESULTS AND DISCUSSION 

Heat- Transfer Rate Distributions Over Cones in Heliiun Test Geis 

The measured heating rates in helium test gas are compared with 
theoretical values for both the cone models in figures 12-15 • The 
measured heating rates for the h0° cone at zero angle of attack are 
compared with theory in figure 12. As discussed earlier, the velocity 
obtained from the pressure film data between stations 27 and 31 
(given in table I) was used for the theoretical computations; the 
other frec-stream quantities were the same for both the cases. It is 
apparent from figure 12 that the theoretical computation using velocity 
obtjvined from pressure film data agtaa'c with Uie measured values within 
6 percent. The measured heating rates for 30° cone at zero angle of 
attack shown in figure 13 also agree with theory (within 7 to 8 percent) 
obtjiined by using U vi-.-ji. However, in this case the measured heating 
rates are generally lower than the theoretical values. This may be 
duo to the fact that the measured heating rates in this case were time 
averaged in a time interval which starts about 55 ysecs from the 


beginning of the I'un, which may have affected the accuracy of 
numerical integration (as mentioned earlier). In figure 1^^, the 
heating rates measured over the 30^ cone at zero angle of attack (a 
repeated run) seem to agree with theoretical values corresponding to 
a velocity (U^) 10 percent more than the one deduced from pressure 
film data (table I). The measured heating rates on the leeward side 
of 30^^ cone at +8*^ angle of attack are compared with theory in figure 15. 

The agreement is quite good in the early portion of the cone stirface, but 
in the aft portion, the agreement is within 10 percent. This is probably due 
to a iittlo more noise in the measured surface temperatures (toward the aft 
portion of tlie cone) than normally existed in other runs. Considering the 
uncertainty in the measured flow velocity, the agreement of the measured 
heating rates over both the cones with theory is considered satisfactory. 

Heat-Transfer Rate Distributions Over Cones in Air Test Gas 

The measured heating rates in air test gas are con5)ared with 
theoretical values for both the cone models in figures l6-20. In figure 
l6, the measured heating rates over the two cones at +12^ angle of 
attack (with gages on leeward side) are given. The theoretical value 
for ^40^ cones are about 5 to 6 percent more than the measured values 
in the beginning and 6 to 8 percent less toward the end of the cone. 

However, no theoretical values were available for 30^ cones. In figures 17 
nnd l8, the measui*ed heating rates over 30 and cones at zero angle of 
attack are compared with theory. The measured stagnation-point heating 
rate using a spherical blunt-nose model (0.025l4~m nose radius) in a 
repeated run is also compared with theory in figures 17 and l8. The 
measured stagnation-point heating rates agree with theory within 2 
percent. However, the measui'ed heating rates over the cone surface 
for both the models are consistently higher than the computed values, 
the differences being 5 to 15 percent in the beginning and increase to 
as much 80 percent toward the end of the cone surface. In figures 
19 <’uid the moa:uired heating rates at angle of attack with gages 
on the windward side, are consistently much more (between 15 to 60 
percent) than tlie computed values. Tlie measurements, which are believed 
to be accurate to about 5 percent (appendix D) thus disagree with the 
oomputfition by too much to attribute to measurement accuracy. At 
present no leasonable explanation (?xists for this discrepancy, it is 
observed from fig\ir<-'s 16-20, that the discrepancy between t>icoiy and 
measurements inezvases as the effective cone oizgle increases. 


CONCLUDING RKMARKS 

Heating-rate measurements over spherically blunted cones have been 
mvade in the expansion tube facility using air and helium test gases. 

The agreement between theory and mcasui’cmonts in the case of helium test 
gas is satisfactoiy . In the case of air test gas, the measured stagnation- 
point heating rate agrees with theory within 2 percent, but the measured 
heating rates over the cone surface on both the cone (liO and 30 ) models 
arc consistently higher than the theoretical values. Tlie discrepancy 
between calculation and experiment increased as the local surface angle 
of attack increased. Tlie remion for the discrepancy is not presently known. 
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APPENDIX A 


CALIBRATION OF HEAT-TRANSFER GAGES 

1/2 

The quartz substrate's physical properties parameter g[g = (kpCp) ' ] 
was determined by measiu:*ing the stagnation -point heat-transfer rate 
to a 0.025^-m radixis spherical model moianted inside a calibration shock 
tube. A thin-film (palladiim) heat-transfer gage was fabricated near 
the stagnation point using the same method as the one used for the heating- 
rate measurements over blunt cones in the expansion tube facility. The 
temperature coefficient of resistance for this gage was measured using the 
same procedure as was used for other gages. The value of a for this gage 
was O.OO2U612 n/R K. 

The calibration shock tube consists of a driver tube (O.Pl^ui long 
and 0.078 ralD) and a driven tube (l.829-m long and 0.15^ mlD) Joined by 
a 0.305~ni long transition piece at the diaphragm station. A polyester 
film diaphragm of 20-ym thickness was used. A shock wave of - 1.5 

was created by using helium as the driver gas and air at ’ 

P2 = 160 torr as the test gas in the driven tube. The shock speed was 
evaluated by measiuring the time of arrival of the shock wave given by a 
counter reading which was started and stopped by shock tube wall pressures 
sensed by two pressure transducers mounted at 0.5336 m apart towards the 
end of the driven tube. The wall pressures which were ’jsed to trigger 
the counter were also monitored in an oscilloscope. The scopes were 
calibrated for gain factor as well as for sweep rate. The pressxure 
transducers were calibrated using a pressure pulse generator. The shock 
speed was also deduced from the measured pressure ratio (pg/pi) as well 
as from the pressure film data. The accuracy of each method is discussed 
later in this appendix. The output from the stagnation-point heat-transfer 
gage was fed into an analog network ( fippendix b) and the output from the 
analog was measui’ed on an oscilloscope. The gain factor for the analog 
was set at 100 and the initial voltage across the gage. Eg, was set at 
1 volt. 

The theoretical expression for the stagnation -point heat-transfer 
rate is given by reference I6; 


= 0.763P;°-® 


R dU 
n e 

Ug dS 


1/2 


(1-A) 


where the Prandtl number Pj. = 0.71 for air at low temperatures. Since 
the flow velocity U2 behind the incident shock wave is subsonic 
in the present case an appropriate nondiraensional velocity gradient 


R dU 

— -rs— given in reference 17 has to be used. The subsonic flow velocity 

Ug is the free-stream velocity for the spherical model. Equation (l-A) 
c^ be expressed as a function of flow Mach number = Ug/a^) and 

(p^/R )l/2 where p^ is the initial pressure in the driven tube and R^ 



APPENDIX A 


is the model nose radius. The variation of the parameter 
flow* Mach number M2 is shown in figure 1-A. 

The stagnation -point heating rate was measured at an initial pressure 
of 160 torr in the driven tube and the incident shock Mach nmber M. 
was around 1.55* The same condition was repeated several times and an 
average value of g was deduced from these repeated runs. The procedure 
for deducing g is illustrated here for a single run and all the data 
obtained in the repeated runs are tabulated in this section. Run No. 11 , 
p^ = 160 torr; R = 0.0254 ra; analog output measured from scope: 10.29 

rav; room tenopera?ure: 25^0; speed of sound a^^ = 346.37 m/sec. Time 
interval measured by counter: 98I psecs * Ug ^ ” 543*9 m/sec. 1 ~ ^•57> 
M2 = 0.66 from charts. ’ 

From figure l.A: 


*^st ^ watts/m^. 


(2-A) 


From the analog heat-transfer relation ( appendix B) 




= 1.975 ^ 10 ^ for tlie analog circuiti 


E = 1000 mv 
0 


D = 100; u = 0.0024612 fl/Q K, 


with the measured value of V = 10.29 mv. 


q ^ = I.OU2 X 10“^(6/a) 
3t 


(3-A) 


Equating (2-A) and (3-A) gives 


IT 


. -s> 


^ i - , ‘ ^ -1 - ' 
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( 0 /o) = 5.7807 X 10^ and 3 «* lk 22 . 7 k w/t/m^OK. 


The difference from the assumed value of 1512.5 is -5.93ltO percent. 

The sh. ck speed Ug 2. was also deduced from the wall-pressure film 
data and gives a shock Mach number Ms^i of 1.575, and the co rrespon ding 
flow Mach number Mg is 0.665 from charts. The value of q . j/il /p. from 
figure 1-A is 2.1*5 x 10^. Then the 8 value deduced from * 
the equations already given is ll*52.3 which is 3.97 percent less than the 
assumed value of 1512.5. 

Since the pressxire pg behind the shock wave was also measured, 
the shock Mach nxamber Mg^l was also inferred from the pressure ratio 
P2/P1 from charts. In this case, pressure behind the shock wave wo’ild 
be pg = 5.56 + 3.091*7 = 6.651*7 psi and p^ = 3.095 psi. Hence, P2/P1 * 
2.797. I ^en M g = ^*59. from charts. The corresponding is O.60 
and qg^rRjj/Pj^ =*2.6 x lo*^ from figure 1-A. Then the 8 value would 

be 151*1.27 which is 1.9 percent more than the assumed value of 1512.5. 

The 8 values deduced by using the preceding method for five repeated 
runs are tabulated below. 


Pjl^ = 160 torr, = 1000 mv, D * 100, a » 0.002l*6l2 ft/0 K 


RUN 

NUMBER 

«3.1 

«2 


9st 

(watts /m2) 

(mv) 

(B/o) 

w/t/m% 

% diff 


1 * 

1.568 

0.660 

2.1*xlol* 

6.023x1oI* 

f 

5.91*7x10^ 

11*63.7 

-3.2 

1 

2 

1.582 

0.672 

2.52xloJ; 

6.32 xior 

10.0 

6.21* xlOc 

1535.5 

+I.5I* 


3 

1.585 

0.672 

2.52x10 

6.32 xlO 

+ 

6.21* xlo’ 

151*5.8 

+I.5I* 


1 

1.578 

0.665 

2.1*5xloi* 

6. 11*9x10 j| 


5.61*3x10^ 

1388.8 

-8.2 

3 

2 

1.602 

0.690 

2. 70x10 J 

6.776x10, 

10,76 

6.219x105 

1530.6 

+1.2 


3 

1.595 

0.685 

2.65x10 

6.651x10** 


6.10j*xl0^ 

1502.3 

-0.67 


1 

1.565 

0.655 

2.35xloi| 

5. 89 8x10 J* 

+ 

5.1t30xlo5 

1335.3 

-11.7 

9 

2 

1.575 

0.665 

2. 1*5x10 J 

6.ll*9xlof 

10.7l» 

5.656x105 

1392.1 

-7.96 


3 

1.610 

0.695 

2.80x10 

7.028x10^ 


6.1*6 xlo5 

1591.1 

+5.2 


1 

I.57U 

0.662 

2. 1*50x10 

6.ll*9xlo{* 

t 

5.50 x 1 q 5 

1353.7 

-10.5 

10 

2 

1.583 

0.672 

2. 520x10 J 

6.325x107 

11. Ol* 

5.567x105 

1392. U 

-7.9 


3 

1.605 

0.690 

2.700x10 

6.776x10** 

+ 

6.06lxl05 

11*91.7 

-l.lj 


1 

1.570 

0.660 

2. 1*00x10 j» 

6.023xloj* 

t 

5. 778x105 

1U22.2 

-5.97 

11 

2 

1.575 

0.665 

2.1*50xl0f 

6. 11*9x10 J 

10.29 

5.901x105 

11*52.3 

-3.97 


3 

1.590 

0.680 

2.600x10** 

6.525x10 

1 

6.262x105 

151*1.3 

+1.90 


*1; M . deduced from counter reading 
s 

2: M , deduced from wall pressure film 

S 9 J. 

3: M , deduced from measured shock pressure ratio 

S 9 X 

Average value of 8 from five runs deduced from shock pressure ratio (item 3 
in the table) is 1532.1* W/T/m^K, which is 1.3 percent more than the assumed 
value of 1512.5. Standard deviation from run to run in the case of 3 is 
2.5 percent from the mean value. 
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From the three preceding methods of obtaining the shock speed, it 
is apparent that b/ using the measured pressiire ratio across the shock 
wave would result in better agreement with the value of g already used 
by other investigators. In order to eliminate the possibility that this 
may be an isolated single cose, the same rm weis repeated 5 times and 
the same trend exists consistently in all the repeated runs as shown in 
the preceding table. The shock speed deduced by counter reading is 
probably affected by two factors: the trigger sensitivity of the counter 

(may not be exactly same for both start and stop gates) euid the initial 
speed of sound a^^ which was computed by using the outside room 
temperatiu:e which may not be the same as the gas temperature inside. For 
the case of shock velocity deduced from the wall pressure film data, the 
probable error from the counter is avoided, but the error due to incorrect 
initial temperature still affects the deduced shock Mach number, in 
the case of shock velocity deduced from the measured pressxire ratio 
(p 2 /pi), both errors are avoided. In this case, we depend on only the 
initial pressure pi which is measrired accurately using an absolute 
vacuum gage and the pressure transducer was calibrated using a pressvire- 
pulse generator. 
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APPENDIX B 


ANALOG NETWORK 

It was realized by Skinner (ref. 20) and Meyer (ref. 21) that electrical 
analogs of a heat-conducting slab could be constructed, which if driven by 
a voltage proportional to surface temperatiare , would yield a signal 
proportional to the heat-transfer rate. Some modifications to these basic 
circuits have been proposed in reference 13. The present analog circuit 
is essentially based on the modifications proposed in reference 13. 


Analogy Between Heat and Current Conduction 

The basis for constructing an electrical analog of a heat-conducting 
slab arises from the similar forms of the equations for heat conduction and 
diffusion of electrical charge in a system in which there is uniformly 
distributed resistance and capacitance. Because of the difficulty of 
constructing an electrical model with continuous variation of capacitance 
and resistance, it is usual to approximate this by employing a network 
of lumped constants. The heat-transfer gage substrate, which is assumed 
to be semi-infiniLe, may be conoi dered to be consisting-'; of a number of small 
strips of slabs of backing material. Each slab is considered as a 
capacitor and resistor. Surface-temperature change corresponds to the 
applied voltage to the network, the corresponding change of current in 
the network corresponds to the heat-transfer rate. The commonly used 
T-section network and the analof"/ botwe<.'n heat and current 
conduction are shown in figure 1-B. The current through the network is 
measured as a voltage drop across the first resistance in the network. 

Then u:>*ng the relation between the heat transfer and the current (fig. 1-P), 
an oxpreu'nion for hoat-t rann for rate can be derived as 

2V (t) fkPc 1 

^ od" [_Y(C J 

O V- 


For a given q(t) the analog output voltage Va(t) is proportional to 
(RC)T/2^ whereas the response time of the analog circuit is given by RC 
so that, in general, it is necessary to make a compromise between fast 
response and adequate output signal level. The network must be made 
long enough so that over the time during which measurements are being 
made, the system behaves as though it were semi-infinite. 

Circuit diagram of analog network and the preamplifier is shown in 
figure 2-B. It consists of 30 RC units with 10 kfi and 390 pf os 
resistance and capacitance in each RC unit, respectively. The first five 
RC units are made of compcaients which are of 5-percent accuracy. This 
will ensure accurate voltage measurement across the first resistance. 

With R and C values used, the response time is 3.9 psecs and the analysis 
of reference 13 shows less than 0.1-percent error due to finite length 
of network. 


21 



APPENDIX B 


The prean 5 >lifier in figure 2-B c«isists of a constant current 
source and a voltage amplifier with a maximum gain of 100. Ihe signal, 
amplified through the amplifier, is directly fed to the analog network, 
but the same signal is attenuated by 10 times and measured as T-out 
which is the surface temperature response of the gage. The gain adjust 
is a precision 10-tum potentiometer which can be used to set the gain 
at any value from 1 to 100. The gain corresponds exactly to the value 
set on the gain adjust potentiometer. The current adjust potentiometer 
is used to set any desired cunrent level through the thin film and two 
binding posts are provided to measure the initial voltage across the 
thin film. The maximum limit of voltage change across the thin film 
should not be more than about 100 to 150 millivolts during the test 
region of any run to avoid saturation in the amplifier. This can be 
maintained by properly setting uhe current level through the gage. 




Figure IB.- Analogy between heat and current conduction. 



icure 2B.- Circuit diagram of analog network. 
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ELECTRIC CIRCUIT USED TO MEASURE SURFACE -TEMPERATURE HISTORY 
The electric circuit used is shown in the sketch given below. 



Let (R^ ♦ Rg) = R^. 

Assumptions : 

1. Rj and Rjj » R^ » R 

2. Temperature coefficient of resistance of thin film a is con-:itant. 
The change in resistance of the thin film at any instant dxiring a run« AR, 
is given by 


AR » (B - R^} « R^o(T - T^) « R^oAT 
or 

AT » AR/(oR ) (1>C) 

o 


where quantities with subscript 0 are initial values at t ■ 0; (AR/R^) 
can be determined from the given quantities R^, R^, and the 
measured change in voltage across the thin film during a run. From the 
circuit shown above, the following relations hold; 


AE - {E - E^) 
o 


(2-C) 


«o"«b 




(3-0 
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E = IR 



(l*-C) 


where E^^ is the battery voltage. 

From equations (2-C), (3-C),and (U-C) we have 



(5-C) 


Using relation given in equation (3-C), equation (;>-C) can be written as: 


AE - M. 

R “ E 
o o 


^ - ^o> 

|_1 - (Ey(E^ - E„))^ 


(6-C) 


Using equation (l-C), the change in surface temperature is given by. 


^ 1 

o J 

AE is the measured change in voltage across the gage. E and E^ i- 
the initial voltage across the gage and the batter;* voltage, respectively. 
The relation given in equation (T-C) was derived by Wilfred Friesen of 
Langley Research Center. 
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ERROR ESTIMATION IN HEAT-TRANSFER MEASUREMENT 

To estimate the absolute error involved in measuring the heating 
rates, the principles outlined in reference 26 are used. Some of the 
important principles are noted below. 

(a) All the errors are classified either independent or dependent. 

(b) In the case of Independent errors, the logical way of adding 
the separate contributions is to take the square root of the sum of 
their squaies. This procedure is considered too conservative and some 
research workers reduce the total error, perhaps by one third. 

(c) In the case of dependent errors, the rule is to replace the 
square root of the sum of squares addition by algebraic addition. 

Uie heat transfers tnat were measured by using electric analog network 
(appendix C) were deduced by using the following expression 


q = 


E oD 
o 



(1-D) 


The error due to the finite length of the analog network is about 0.1 
percent ( appendix B) and the value of RC can be taken as quite accurate 
since the values of R and C in the network are of 0.1-percent accuracy. 
Then, the change in the heating rate can be written from equation (l-D) as 


AV 

Aq \ + M. M 

q ■ 3 “ a 


The error in meeisuring by using the recording system is estimated 
to be less than 2 to 3 percent. The error in 3 can be taken as 2.^ 
percent ( appendix A) . The value of o of the thin film was obtained 
by normalizing with the resistance of the gage measured at ambient 
temperature dxiring calibration. But the initial temperature of the 
thin film Just before the run is made will be different since the 
ambient temperature will vary and in addition, the thin film gets heated 
due to constant current passing throiigh the film. The error introduced 
due to these effects is taken conservatively as 3 percent. 

Then 
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2 

+ 2.5 + 


*<•6 pei'oont 


TJuir., a conservative error estimate would be around 5 percent. 

to deduced iVom numericjil integration are sublected 

to the s tune errors os noted above, but on additional error in selecting 
hi time intep-ation steps is introduced. This additional error is shown 
Ihi^ronortf "" selecting different time steps os noted elsewhere in 


f 
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TABLE I.- MEASURED FREE-STREAM VELOCITIES 


Test gas: Air 


Run 

number 

Velocity at tube exit 

Velocity between stations ^ 

and @ (1.551* m) 

(from st.-line fit to 
the data) 

From coxHiter 
reading* 

From HT film 
data 

From pressure film 
data 


(m/s) 

(m/s) 

(m/s) 

(m/s) 

EC IU77 

5061 

5090 

5096 

5062 

EC 1U78 

5092 

5182 

5062 

5096 

EC ll»79 

5130 

5212 

5215 

5129 

EC ll*9** 

51li* 

511»6 

5215 

5096 

EC 11»95 

5176 

5198 

5251 

5096 

EC 1U96 

5157 

5251 

5378 

5129 

EC IU97 

5203 

5323 

531*1 

5215 

EC IU98 

5026 

511*6 

- 

- 

Test f:as : 

Helium 

^Cotmters are 

triggered by heat-transfer gages. 

EC lk73 

61^10 

6396 

6396 

6787 

EC IU75 

63i*2 

621*2 

6318 

681*7 

EC 1U76 

6376 

6267 

6292 

6728 

EC IU91 

61*7J* 

6318 

63U5 

6787 

EC ll»99 

6253 

6252 

- 

- 

EC 1500 

61*15 

6267 

6585 

6670 


Station .ocations: (2^ : 19.732 m; : 21.286 m; tube exit: 21.610 in. 

Mid-point location between St. and ^ : 1.1 m from tube exit. 
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TABLE II.- TYPICAL o VALUES FOR THIN- FILM GAGES 


Gage 

numbei’ 

Gage resistance 
at room tenq?. 
( 2 k . k ^ C ) 


0 values (n/£J°K) 


38.39°C 

52.2°C 

66.2°C 

79.5°C 

1 

- 


- 

- 

- 

2 

77.8(0) 

0.002364 

0.002373 

0.002377 

0.002367 


76.9 

0.002366 

0.002373 

0.002379 

0.002368 

k 

76. U 

0.002361 

0.002367 

0.002373 

0.002363 

5 

77.7 

0.002361 

0.002366 

0.002373 

0.002363 

6 

78.8 

0.002362 

0.002366 

0.002374 

0.002364 

7 

80.5 

0.002363 

0.002369 

0.002377 

0.002367 

8 

81.9 

0.002358 

0.002367 

0.002375 

0.002366 

9 

83.7 

0.002362 

0.002370 

0.002379 

0.002370 

10 

9 h .6 

0.002330 

0.002310 

0.002323 

0.002315 

11 

84.8 

0.002359 

0.002370 

0.002380 

0,002370 

12 

83.2 

0.002362 

0.002372 

0.002383 

i 

0.002373 

13 1 

83.5 

0.002365 

0.002370 { 

0.002378 

0,002369 

ll» 

83.6 

0.002367 

0.002369 

0.002379 

0.002370 

15 

82.6 

0.002365 

0.002366 

0.002377 

0.002372 

l6 

87.6 

0.002351 

0.002372 

0.002380 

0.002371 

Standard deviation 
(from the mean value) 

OM 

0.65!? 

0.6!? 

0.6!? 
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TABLE III.- INITIAL TEST CONDITIONS 
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TABLE IV.- MEASURED FREE-STREAM CONDITIONS 


Test -as : Air 


Run 

number 

Pt o 
(N/m^) 

Poop 

(N/m^) 

'^w 

(•’K) 

9 st* 

(W/m^) 

U27-31 

(m/sec) 

EC 11+77 

108058 

I68I+ 

325 

— 

5062 

EC i 1+78 

11371*7 

1620 

319 

— 

5096 

EC ll +79 

116170 

1606 

356 

— 

5129 

EC ll* 9 l* 

113800 

1682 

31+5 


5096 

EC ll *95 

111031 

1722 

359 


5096 

EC i 1+96 

119102 

1776 

360 

— 

5129 

EC lh 9 T 

— 

2052 

380 

1.31+8x10'^ 

5215 

EC 11+98 

— 

2000 

380 

1.277x10’^ 

5026 

Test gas 

: Helium 





EC 11+73 

52299 j 

958 

3 ll* 

— 

6787 

EC-ll +75 

51965 

863 

328 

— 

681+7 

EC 11+76 

1+9678 

807 

313 

— 

6728 

EC 11+91 

51299 

831+ 

317 

— 

6787 

nc 11+99 


815 

333 

7 . 69 x 10 ^ 

— 

EC 1500 

— 

963 

330 

8.3x1o'^ 

6670 


«Nose radius Rj^ = 0 . 025 ^ m for stagnation-point 
heating -rate measurements. R^ = 0.00635 m for 
blunt-nosed cone models (as shown in fig. 5). 
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TABLE V.- COMPUTED FREE-STREAM FLOW CONDITIONS (REF. ll) 


Test gas: Air 


Run 

number 

Density 

(kg/ra3) 

Temperature 

(°K) 

Enthalpy 

(in^/sec^) 

Compressi- 

bility 

Specific 
heat ratio 

Mach 

number 

Reynolds 
number 
per meter 

EC IU 77 

4.35x10"^ 

1349 

1 . 456 x 10 ^ 

1.010 

1.313 

7.10 

4.38x10^ 

EC 11*78 

4.52x10”^ 

1248 

1 . 336 x 10 ^ 

1.010 

1.319 

7.41 

4.81x10^ 

EC ll*79 

4 . 56 x 10 "^ 

1227 

1.311x10^ 

1.010 

1.320 

7.52 

4.94x10^ 

EC ll*9l* 

4.52x10"^ 

1297 

1.394x10 

1.010 

1.316 

7.28 

4.70x10^ 

EC 1 I 95 

4.41x10"^ 

1361 

1.470x10^ 

1.010 

1.312 

7.12 

4.45 10 ^ 

EC 1 L 96 

4,67x10"^ 

1325 

1.428x10^ 

1.010 

1.314 

7.25 

4.82x10^ 

Test gas 

: Helium 







EC ll*73 

1 . 277 x 10 "^ 

360 

1 

1 . 87 x 10 ^ 

1.000 

1.667 

1 

6.08 

3.59x10^ 

EC 1475 

—4 

9.5x10 

331 

1.72x10^ 

1.000 

1.667 

6.39 

3.75x10^ 

EC 1476 

1.235x10*^ 

313 

1.63x10^ 

1.000 

1.667 

6.46 

3 . 80 x 10 ^ 

FC 1491 

3 .254x10"^ 

318 

1.65x10^ 

1.000 

1.667 

6.46 

3.84x10^ 
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Figure 1.- Schematic diagram of expansion tuhe flow sequence 
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Figure 2.~ Pitot-pressure profile at various distances downstream of tube exit 
P- S 3 .U 5 kN/m^; P^_ * I 6 .OO N/m^; t = 200 \is. 
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Figure 3.- Variation of shock velocity, U along acceleration section. 
Test gas: Helium, Run lUpi. 
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Figure J*.- Variation of shock velocity, U , along acceleration section. 
Test gas: Air, Run ll*'/9. 
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9*- Conparison of hea“ir.g rates obtained by numerical integration of 

shock-tube wall-surface temperature with results of analog circuit. 




Figure 10.- Comparison of heating rates obtained by numerical integration of 

stagnation-point temperature history with results of analog output 








Figure 11.- Continued. 
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i^rure lU.- Cocparison of a«iasured heating rates with theory for helium test gas 
F,un'lU91. 
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'ure _o.- Cor-r.ariccn of .T.v'i.-ured hc-uting rates vith theory for 30° and 1*0° cone 
ciodel.- at +1?° angle of attack. Air test gas. Runs 11^77 and ll73. 







^•asurwl stagnation* point lieating rate 


O 









rates with theory for 30® cone mod 






